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Abstract 

This paper develops the asymptotic matching of an unsteady compressible boundary 
layer to an inviscid flow. Of particular importance is the velocity injection or transpiration 
boundary condition derived by this theory. It is found that in general the transpiration will 
contain a slope of the displacement thickness and a time derivative of a density integral. 
The conditions under which the second term may be neglected, and its consistency with 
the established results of interacting boundary layer are discussed. 


Introduction 

A variety of viscous/inviscid interaction schemes have been developed for a steady or un- 
steady boundary layer coupled with steady or unsteady inviscid flow. For various reasons (e.g. 
expediency, robustness, efficiency), engineering codes have frequently combined unsteady invis- 
cid and steady boundary layer flow solvers. In such cases the coupling typically is treated as 
quasi-steady. The attempt will be made here to show that under the correct circumstances 
this approach is approximately valid. The larger question is the proper way to derive the invis- 
cid velocity boundary condition imposed by an unsteady boundary layer. Heuristic treatments 
appeal to the similarity between the displacement thickness slope derived from boundary layer 
theory and the surface slope derived from thin airfoil theory. Based on this similarity, a time 
derivative of the displacement thickness is sometimes added to the viscous/inviscid coupling. 
This formulation has shown up from time to time, as for example a viscous-inviscid interaction 
model for cascades, [1] or as a recent publication apparently illustrates. [2] In contrast to this 
approach, the attempt will be made here to show that the method of matched asymptotic ex- 
pansions applied to an unsteady interacting boundary layer (IBL) yields only a time derivative 
due to the density variation through the boundary layer. This result was first presented by 
LeBalleur, [3] and later by Bartels [4] [5] [6], and Epureanu [7]. 

The paper will derive the transpiration boundary condition for an unsteady compressible 
laminar boundary layer. The same form will be shown to apply to an unsteady turbulent 
boundary layer. The transpiration boundary condition will be derived for a surface with time 
and spatially varying height using the Prandtl transposition. Reduced forms of the transpiration 
condition will be shown in the limits of incompressibility and steady state, which are then shown 
to merge with other well known interacting boundary layer results. Finally, its relation to the 
IBL theory of boundary layer separation and other IBL methods will be discussed. 
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Laminar boundary layer matching 

A laminar two-dimensional boundary layer is composed of two regions, an inner viscous and 
an outer inviscid but rotational region. The leading order boundary layer equations, matching 
conditions and the resulting injection velocity will be derived using the method of matched 
asymptotic expansions, developed for fluid dynamics by Van Dyke. [8] 

The Navier-Stokes equations are non-dimensionalized with free stream density, velocity and 
characteristic length l that is of 0(1). The length scale on which the present interaction 
is derived is x ~ 0(1). The Reynolds number is defined such that Re = pu^l/p. The non- 
dimensional inviscid flow field density, x and y velocities, pressure and temperature are expanded 
in the following way. 


Pi(x,y,t) + ep 2 (x,y,t) + ■■■ 

(1) 

ui(x,y,t) + eu 2 (x,y,t) H 

(2) 

vi (x,y,t) + ev 2 (x,y,t) H 

(3) 

Pi(x,y,t) + ep 2 (x,y,t) + • • • 

(4) 

T\{x,y,t) + eT 2 (x,y,t ) H 



The scale parameter e for a laminar boundary layer is Re 1,/2 . The vertical scale of the inviscid 
flow is y ~ 0(1). The inviscid mass and momentum continuity equations at leading order are 


Pit + (Pi u i)x + (PlVl)y = 0 (5) 

{piU 1 )t + {p l u 2 1 ) x -\-{p 1 viu 1 ) y -\-p lx = 0 (6) 

(PiVi)t + (Pl«Rl)x + (pivl)y + Piy = o (7) 

Athough the energy equation can be similarly expanded, it is not necessary for the present 

purpose and will not be included. The boundary layer expansions have the following form 

p ~ Ri(x, Y,t) + eR 2 (x, Y, t) + • • • (8) 

u ~ Ui(x, Y,t) + eU 2 (x,Y,t) + • • • (9) 

v ~ eV 1 (x,Y,t) + e 2 V 2 (x,Y,t) + - ■ ■ (10) 

P ~ pi(x,Y,t) + ep 2 (x,Y,t)-\ (11) 

T ~ Q 1 (x,Y,t) + eQ 2 (x,Y,t) + --- 


where the boundary layer scaled coordinate Y = ye 1 is used. The leading order boundary 
layer equations are 


Rit + (R\Ui) x + (-RiVi)y — 0 (12) 

(RiUi) t + (RiUi) x + (Ri\iUi)y +Pi x — (pUiy)y (13) 

Piy = 0 (14) 

The method of matched asymptotic expansions requires that the Taylor series of the inviscid 
flow quantities as y — > 0 is matched with the boundary layer expansion as Y —> oc. [8] The 
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Taylor series of the inviscid expansion is 


P ~ Pi(x,0,t)+yp ly (x,0,t) + ep 2 (x,0,t) + --- 
~ Pi{x,0,t) + e(Yp ly (x,0,t) + p 2 (x,0,t)) + ■■■ 
u ~ Ui(x, 0, t) + yui y (x, 0, t) + eu 2 (x, 0, t) + ■ ■ ■ 

~ ui(x, 0, t) + e (Yui y (x, 0, t) + u 2 (x, 0, t)) + • ■ • 
v ~ vi(x,0,t) + yv ly (x,0,t) + ev 2 (x,0,t) + ■ ■ ■ 

~ vi(x, 0, t) + e(Tni J/ (a;,0,t) +n 2 (x,0, t)) + • • • 
etc... 

Matching of each quantity at successive orders as y — > 0 with the corresponding boundary layer 
quantities as Y —> oo yields the matching at leading order 

Pi_(x, 0,t) = R]_(x, Y, t) 
ui(x,0,t) = Ui(x,Y,t) 

Vi(x,0 ,t) = 0 

etc... 

and 

p 2 (x,0,t) + Yp ly (x,0,t) = R 2 (x,Y,t ) 
u 2 (x, 0,t) + Yui y (x,0,t) = U 2 (x, Y, t) 
v 2 (x,0,t) + Yv ly (x,0,t) = Vi(x,Y,t) 

etc... 

Note that 

Riy 0 ( 21 ) 

as Y — > oo for density to remain finite in the limit. In view of equations 15, 16 and 21, the 
correspondence can be made between the continuity equations in the dual limits as y — > 0 and 
Y — > oo that 


(18) 

(19) 

( 20 ) 


(15) 

(16) 
(17) 


—R.\V\y = R\t + {RiUi) x in the Inner Limit (22) 

= Pit + (Pi u i)x = — (Pi v i)y in the Outer Limit (23) 

By making use of equations 15 and 17, equations 22 and 23 can be written 


RlV 1Y = PlVly 


or finally that 


(24) 


V 1Y (x, Y, t) = v ls (x,0,t) (25) 

as y — > oo. The expressions for velocity matching at leading order, equations 16 and 20 can be 
written 


ui(x,0,t) = Ui(x,Y,t ) 

v 2 (x,0,t) = Vi(x,Y,t) - YV 1Y (x,Y,t) 
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(26) 

(27) 



This is identical to the leading order velocity matching for a boundary layer derived by Van 
Dyke. [9] These two matching conditions are finite in the limit Y — > oo and provide, along with 
expressions for density, pressure and temperature, the necessary inviscid/ viscous matching for a 
compressible unsteady boundary layer. The first of these states the widely recognized matching 
of the rr-component of the inviscid and viscous velocities at the boundary layer edge. The second 
is the boundary layer contribution to the injection or transpiration velocity into the inviscid flow. 
The injection velocity can easily be rewritten into a more usable form. Integrating the boundary 
layer continuity equation 12 yields the result 


p 1 (x,0,t)V 1 (a;,V,t) = i? 1 (a;,V,t)V 1 (a;,V,t) = - [ (R lt + (RiUj^dY (28) 

Jo 


in the limit as Y — ^ oo . This integral is infinite in the limit as Re — > oo. However, by using 
the velocity matching of equation 27 combined with equations 5, 28 and 25, the finite condition 


v 2 (x,0 ,t) 


V l {x,YR) - YV iy {x,YR) 



{RiU x ) x )dY 



1 

Pi 



(. PiUi - R\U\)dY + 



Rit)dY^j 

R^dY) 


is obtained. Using the definition of displacement thickness 6* 



and defining the new term, density thickness 8r 

roc JD 

8r = (1 - — )dY 

Jo Pi 


the resulting injection velocity of an unsteady compressible boundary layer into the inviscid flow 
is 


V2 (x,0,t) 


± ( d{PiU\8*) d(p l 5 R ) \ 
p l \ dx dt ) 


(29) 


This important result presents an asymptotically consistent treatment of the influence of an 
unsteady boundary layer on the outer inviscid flow. It is identical in form to the matching 
derived by Le Balleur [3] and later by Bartels [4], [6], [5] and Epureanu. [7] It states that the 
inviscid flow sees an injection velocity due to the sum of the slope of the displacement thickness 
and the time variation of the integral of density through the boundary layer. 

Now several simplifications of this theory can be found. If the flow is unsteady but incom- 
pressible 

djpM = n 

dt 

and the injection velocity can be written 


v 2 (x,0,t) 


dimS*) 

dx 


(30) 
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where the incompressible displacement thickness is 


poo TT 

6* = (1 - — )dY . (31) 

Jo u i 

This is the widely used incompressible boundary layer velocity matching condition. If the 
density variation through the boundary layer is small or the time scale of the variation is long 
(e.g. low frequency or quasi-steady) such that 


d(PiM 

dt 


~ 0(e) 


the injection velocity can be written 


v 2 (x,0,t) fa 


1 d(p 1 u 1 S*) 
p x dx 


(32) 


This is the widely used quasi-steady interaction. It is a reasonable approximation for many 
subsonic and transonic adiabatic flows. It is identical to the velocity condition for a steady 
compressible subsonic or supersonic boundary layer used by Davis [10]. 


Turbulent boundary layer matching 


In contrast to the two layer structure of a laminar boundary layer, a turbulent boundary 
layer has three layers in the limit of large Reynolds number. Nevertheless, if the necessary 
equations of turbulent theory are identical to the corresponding equations of laminar theory it 
is possible to utilize the same matching conditions between the appropriate regions for laminar 
and turbulent boundary layers. This will be shown here. The present formulation follows 
that of Mellor. [11] (See also ref. [12])The turbulent boundary layer is expanded in the small 
parameter e = ut/uo where u t is friction velocity or turbulent velocity at some designated point 
and uq is a characteristic flow velocity, e.g. free stream mean flow velocity. The outer inviscid 
region is similar in nature to the inviscid region of a laminar boundary layer. It is expanded in 
a form identical to that of equations 1-4. 


p 

~ Pi(x,y,t) + ep 2 (x,y,t) + ■ ■ ■ 

(33) 

u 

~ ui(x,y,t) + eu 2 (x,y,t) 4 

(34) 

V 

~ v 1 (x,y,t) + ev 2 (x,y,t)-\ 

(35) 

p 

~ Pi(x,y,t) + ep 2 (x,y,t) + - ■ ■ 

(36) 

T 

~ Ti{x,y,t) + eT 2 (x,y,t) 4 

(37) 


This expansion follows that of Mellor [11], except that temperature and density are also ex- 
panded. The turbulent stress expanded in ref. [11] in the inviscid, defect and viscous regions 
is not necessary for the present purpose. The first order inviscid equations are identical to 
equations 5-7. 


Pit + (plUl)x + (PlVl)y 

= 0 

(38) 

(PlUl)t + (PlUl)x + (PlVlUl)y + Pl X 

= 0 

(39) 

{PlVl)t + (plUlV^x + (pll’Dy+Ply 

= 0 

(40) 
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The term by term expansions in the defect region follows the form shown in equations 8-11. 


p ~ Ri(x, Y,t) + eR 2 (x,Y,t) + ■ ■ ■ (41) 

u ~ U 1 (x,Y,t) + eU 2 (x,Y,t) + --- (42) 

v ~ eV 1 (x,Y,t) + e 2 V 2 (x,Y,t) + --- (43) 

P ~ pi(x,Y,t) + ep 2 (x,Y,t)-t (44) 

T ~ Q 1 (x,Y,t)+e0 2 (x,Y,t) + --- (45) 

As with the inviscid expansion, this follows that of Mellor [11] with the addition of temperature 
and density expansions. The coordinate normal to the wall is stretched according to Y = ye~ 1 . 
The leading order equations in the defect layer are 

Rit + {RiU\) x + (-RiRl)y = 0 (46) 

(R 1 U 1 ) t + (R l U!) x + (R 1 V 1 U 1 ) Y +p lx = 0 (47) 

Piy = 0 (48) 


The defect layer equations have a zero normal pressure gradient as do the laminar boundary 
layer equations. The defect layer equations are inviscid at leading order and thus, in structure, 
appear as a subset of the laminar equations. In contrast to the inviscid outer region, it contains 
Reynolds stress terms at second order. The viscous layer at leading order also preserves zero 
normal pressure gradient. This is the layer from which the log law for a turbulent boundary 
layer is derived. 

The matching of the inviscid and defect layers is the result of interest. At leading order it 
can be written, following Mellor [11] 


ui(x,0,t) = Ui(x,Y,t ) 

Vi(x,0,t) = 0 

v 2 (x, 0,t) + Yvi y (x,0,t) = Vi (x,Y,t) 

For a compressible flow the matching of density, among other quantities, is also required. At 
leading order 

Pi(x, 0,t) = Ri(x,Y,t) 

The matching conditions at leading order between the inviscid and defect layer for a turbulent 
boundary layer are thus identical to those matching the inviscid with a laminar boundary layer. 
Furthermore, since the mass continuity equations of the inviscid and defect layers for a turbu- 
lent boundary layer (equations 38 and 46) are identical to those for a laminar boundary layer 
(equations 5 and 12), the end result is that the injection velocity due to a turbulent boundary 
layer will have a form identical to that for a laminar boundary layer. 

Boundary layer over an airfoil and the transfer of the 
transpiration boundary condition to the X-axis 

This formulation follows that in Davis and Werle in which the baseline to which the boundary 
layer equations are transfered is the rc-axis. See [13] , [14] and [15] for details. The geometry 
and coordinate systems used in the Prandtl transposition are shown in Figure 1. 
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Figure 1: Geometry for Prandtl transposition 


The Prandtl shift of the boundary layer to a boundary layer along the x-axis requires defin- 
ition of the normal coordinate 

Y = Y + f(x,t) 

along with the normal boundary layer velocity 

V(x, Y , t ) = V 1 {x, Y, t) - f^x, Y, t ) - f t 

The continuity, momentum and energy equations can be shown to have the same form when 
transformed. The matching condition (equation 29) as Y — > oo is written 

v 2 (x, 0, t) = V 1 - YV\ Y = V- YV 9 + uj x + f t - fV 9 


But since it can be shown that 


V - YVy 


JL ( d (Pi u i S *) , 

p 1 \ dx dt ) 


and 

djpiuj) d{pj) 
dx dt 

as Y — > oo, the resulting transpiration boundary condition is 


uifx + ft - fVy = — 

Pi 


v 2 (x,0,t) 


f d{p]Ui{5* + /)) + /)) 

p 1 \ dx dt 


This is the transpiration condition for a surface with height / transformed using the Prandtl 
transposition. For small time variation in the density and small spatial variation in p x u \ this 
can be approximated as 


v 2 (x,0 ,t) 


1 

Pi 


djpjUiY) d(pM 

r\ r 


dx 
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dt 


+ Mi f x + ft 



Consistency with unsteady interactive boundary layer 

theory 


It is useful to show that a proposed theory is consistent with other established results. 
Comparisons with the results of steady boundary layer asymptotics have been made earlier in 
this paper. Additional comparisons can be made with theories developed for unsteady boundary 
layer interaction. The present paper predicts for 8 rt — 0 that the unsteady laminar transpiration 
velocity for a problem with x ~ 0(1) is given by the slope of the displacement thickness only, 
i.e. 


v(x, 0, t) = 


d(u8* 


Re" 1 / 2 dx 

This is consistent with the unsteady laminar IBL on the short length scale of a separation. On 
this scale an incompressible viscous lower deck interacts with the outer inviscid flow through 
the pressure-displacement relation 


p(x,t) 


i r tm* 

7T J - oo X -S 


for a subsonic incompressible or compressible free stream involving the slope of the unknown 
displacement function A(x. t). [16] For a supersonic free stream the pressure-displacement rela- 
tion is p = —A x (x,t). [16], [17], [18]. This is also the interaction of an unsteady boundary layer 
separation at the incipient development of Tollmien-Schlichting waves. [14] Peridier, Smith and 
Walker apply this IBL theory to an x ~ 0(1) problem of a laminar incompressible boundary 
layer interacting with a vortex. [19], [20] To simplify somewhat the form, their equations 72-74 
(ref. [20]) are rewritten here for a flat plate boundary layer in the absence of external excitation. 
The interaction at leading order is 

u(x,0,t) — ui(x,0,t) + • • • = 1 -I —TJ? f — ds (49) 

vr Re 1/2 J- oo x ~ s 


where the outer flow injection is induced by the displacement slope, i.e. 

F{x,t) = ^ («i(x,0,i)<5*(a;,i)) . (50) 

The displacement slope F(x,t) is reproduced from ref. [20]. At the same time it and the 
velocity- displacement relationship (equation 49) are directly obtainable from the results derived 
in the present paper for a flat plate boundary layer in which 5 j> = 0. 

The present theory is in agreement with the viscous-inviscid overlay or "defect formulation" 
of Le Balleur. [3] The velocity boundary condition to the inviscid flow is derived for an unsteady 
compressible flow. In the notation of ref. [3] the transpiration velocity is w(x,0,t), boundary 
layer edge velocity and density are p and u. Using the inviscid and viscous continuity equations, 
the result presented on page 24 of ref. [3] is 

d f°° d 

pw(x,0,t) = — j ( p-p)dz + —J (pu — pu)dz 

where the quantities with the overbar are due to the boundary layer. Epureanu [7] presents 
the transpiration flux Q u 

r , _ d(p e Sp) d(p e u e S*) 

VM — qr 1 x • 

at ox 

8 



Other than notational differences (S p = Sr in the present notation), these expressions are equiv- 
alent and identical to equation 29 of the present paper. 

Concluding remarks 

This paper attempts to clear up misconceptions regarding the proper approach to coupling 
an unsteady boundary layer and inviscid flow solver. The general form of the transpiration 
velocity has been derived for laminar and turbulent boundary layers. Various simplifications of 
the transpiration velocity are also derived. The most general form of transpiration as well as the 
simplified forms are all demonstrated to be consistent with the results of other well established 
methods and are fully consistent with the asymptotic theory of unsteady incompressible and 
compressible boundary layer separation. 
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